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Abstract 

Numerical simulations shed light on control of shear elastic wave propagation in plates structured 
with inertial resonators. The structural element is composed of a heavy core connected to the main 
freestanding plate through tiny ligaments. It is shown that such a configuration exhibits a complete 
band gap in the low frequency regime. As a byproduct, we further describe the asymmetric twisting 
vibration of a single scatterer via modal analysis, dispersion and transmission loss. This might pave 
the way to new functionalities such as focusing and self-collimation in elastic plates. 

PACS numbers: 43.20.+g, 43.35.+d, 77.65.Dq 
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Over the late nineties, structured materials have attracted particular attention due to 
their functionalities regarding waves’ control. Regardless of the wave vector, the direction 
of propagation and polarisation, phononic crystals and locally resonant sonic crystals can 
prohibit the propagation of elastic waves by means of destructive interferences due to Bragg 
scattering and local resonances of each individual scatterer, respectively DEI . Moreover, 
waves can easily be channelled through waveguides j3], slowed down or trapped to enhance 
interaction between photons and phonons for example |1]. Furthermore, locally resonant 
sonic crystals can be viewed as acoustic metamaterials, counterparts of electromagnetic 
metamaterials, which were introduced by Pendry and his colleagues to enhance coupling 
between resonators [5]. 

The hrst studies were dedicated to bulk waves where no supplementary boundary condi¬ 
tion were taken into account. The concepts discussed above have been quickly transposed 
to surface guided waves and Lamb waves mm- In 1998, Tanaka and Tamura calculated 
the band structures to investigate the dispersion of surface elastic waves in a two dimen¬ 
sional semi-inhnite phononic crystal [S]. The directional and complete surface band gaps 
in drilled-hole phononic crystals have been presented by Wu et ah and Benchabane et ah, 
respectively [91110]. The locally resonant mechanism was compared to Bragg scattering for 
surface elastic waves a few years later m- 

Unlike for surface elastic waves where very few studies have been conducted, many papers 
were published on waves propagating in structured plates; especially flexural waves where 
both analytical approximations (Kirchhoff, Mindelin) and exact numerical solutions were 
reported. The zero frequency band gaps in arrays of inhnite conducting wires na and in 
pinned plates [T3|, lensing in platonic crystals m or low frequency band gaps in locally 
resonant sonic crystals with stiff pillars [151 Il6] , heavy pillars combined with soft rubber Id 
or Helmholtz-like resonators [T8| are examples of those. 

In the case of plates with resonators, many researchers explored different conhgurations 
to engineer band structures by tailoring the overall shapes of the device. For instance, 
Assouar et al. ED and Wang et al. [22] studied numerically the enlargement of locally 
resonant acoustic band gaps and the propagation of Lamb waves in two-dimensional sonic 
crystals based on a double-sided stubbed plate. Wang et al. implemented the optimized 
design of alternate-hole-defect on a silicon PC slab in a square lattice. Recently, Zhang et 
al. introduced a spiral-like resonator as a unit cell to obtain large band gaps in the low 
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Figure 1: Schematics of three configurations used for investigating the propagation of shear 
phonons in a periodic array of inertial resonators linked to a 2 mm-thick plate through short tiny 
ligaments, a) Both the plate and the resonators are made of duraluminium and the ligaments 
are aligned, b) The central inclusion is made of tungsten (heavy material) and connected to the 
host plate through aligned ligaments, c) Same as b) but now with ligaments tilted by 20 degrees, 
d) Schematics of the periodic structure used to achieve collimation and focussing. The central 
inclusions are made of tungsten and connected with ligaments tilted by 20 degrees / -20 degrees in 
order to facilitate the change of the direction of the wave propagation at the resonance frequency. 

frequency range [23] . 

In this letter, we describe a structure composed of inclusions, heavier than the host plate, 
sustained by very tiny ligaments. We first demonstrate that such a configuration enables 
a strong attenuation of every elastic wave polarisation in a wide frequency range compared 
to its equivalent monolithic counterpart. We present afterwards the asymmetric twisting 
vibration mechanism induced by tilted ligaments and its aptitude to focus elastic shear 
waves. Such structures have been reported previously in the case of bulk waves but were 
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Figure 2: Band diagram of elastic wave propagation within platonic crystals with inertial res¬ 
onators corresponding to the three conhgurations in Fig. 1(a) (upper panel), Fig. 1(b) (middle 
panel) and Fig. 1(c) (lower panel). F, X and M denote the vertices (0,0), (vr/d, 0) and (vr/d, vr/d) 
of the first Brillouin zone, where d is the pitch of the array. 
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limited to the demonstration of band gaps [191120] • 

The nnit cell of the model considered here is a 3.5 mm radins cylindrical inclnsion em¬ 
bedded in a perforated plate, snstained throngh tiny ligaments. Inclnsions, ligaments and 
the host plate are all 2 mm thick. We compare three conhgnrations, namely, the mono¬ 
lithic dnralnmininm plate with aligned ligaments, tnngsten inclnsions in a dnralnmininm 
plate with aligned ligaments and hnally, tnngsten inclusions in a duraluminium plate with 
ligaments tilted by 20 degrees (from top to bottom panels in hgure [^. Since the whole 
structure is of constant thickness, the plate could be easily incorporated in plate devices 
and is a good candidate to replace freestanding membranes. The physical parameters used 
in these simulations are ptun = 17800 kg/m^, Xtun = 166 GPa and ptun = 129 GPa for 
tungsten and pai = 2790 kg/m^, Aaz = 56.2 GPa and paz = 28.7 GPa for aluminium. We 
hrst perform the band structure calculation in the hrst Brillouin zone by imposing Floquet- 
Bloch conditions on opposite sides of the unit cell, see Fig. In hgure [^a), it is shown 
that the presence of the resonators does not disturb the wave propagation up to 60 kHz. 
The frequency range [60 kHz, 100 kHz] is characterized by the presence of coupling between 
the inclusions’ resonance and in-plane coupled pressure and shear waves. The resonance 
frequency observed at 70 kHz of the twisting mode can be estimated by the formula m ■■ 


/ = 7r-< 


N(^X -|- /^)/^ 


N 




( 1 ) 


27r'^ 3(A + 2p)M ^ {k - 0^)3 

Here M = R'^tt p, denotes the mass of the inner resonator and R being the inner 
radius, eh = 0.2 mm the ligaments width, bi — ai = 0.5 mm the ligaments length 
and iV = 4 the number of ligaments involved. When the mass is connected to the plate 
through one ligament, both the calculated and the computed resonance frequencies are equal 
to 37 Hz. It should be noted that in the above formula p has the dimension of kg/m^ since 
in-plane waves do not sense the third direction. Unlike the monolithic plate where both 
the ligaments and the inner mass are made of the same material as the host plate, this 
formula fails to estimate the resonance frequency when different materials are considered. 
For the rest of the document, only results based on hnite element numerical simulations are 
presented. The inner duraluminium resonator in the unit cell has been replaced afterwards 
by a tungsten one, the same geometrical properties hold (hgure [^b). In hgure [^b, we depict 
the dispersion curves of this structure. It is shown that a complete band gap occurs between 
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Figure 3; Illustration of the twisting modes (a- at 30 kHz, b- at 32 kHz ) in the case of aligned liga¬ 
ments and ligaments titled by 20 degrees of tungsten-inclusion-based platonic crystal. It highlights 
the energy transfer from the twisting mode to the plate. The reader may refer to supplementary 
materials for better viewing. 

40 kHz and 60 kHz. Moreover, one can notice the presence of a flat band at 30 kHz. This 
corresponds to a rotational mode, difficnlt to conple with the continnnm becanse of the 
heavy core resonator. In order to facilitate this inertial resonance, we tilt each ligament by 
twenty degrees. In hgnrej^c) we draw the dispersion cnrves for snch a strnctnre. It is worth 
noting that these bands match nearly perfectly the dispersion cnrves obtained in the case 
of aligned ligaments (hgnrej^b)) except near the resonance between 30 kHz and 40 kHz. 
Indeed, a strong conpling between the rotational mode on the one hand, the shear wave 
and the compressional wave on the other hand, is snggested by nnmerons level repnlsions 
between bands. We fnrther illnstrate this hybridization between the continnnm and the 
local resonance of the inclnsion in hgnre It can be easily seen that the resonator in the 
basic cell twists independently withont indncing any collateral vibration of the host plate in 
the case of aligned ligaments. By contrast, the resonator commnnicates more energy to the 
host plate when the ligaments are tilted. The phenomenon is best seen with movies of the 
eigenmodes’ vibration (see snpplementary materials) [2S]. This strong conpling between the 
host plate and the resonators is accompanied by a redirection of shear waves in a privileged 
direction cansed by the intrinsic asymmetric shape of the resonators. This observation can 
be made regardless of the direction of propagation. This makes these kinds of resonators 
good candidates for focnsing, collimation and prohibition of waves, as we shall now explain. 

Fignre l^snmmarizes the aforementioned fnnctionalities. In each snb-hgnre, we plot the 
magnitnde of the held, its phase shift and its iso-snrfaces in order to highlight physical phe- 
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Figure 4: Energy cartography, induced phase shift and iso-surfaces of shear in-plane wave in the 
presence or absence of a seven-column platonic crystal made of rotators with tilted ligaments. It 
underlines different functionalities that may be achieved at different frequencies, a) Endoscope 
effect, b) focusing and c) total mirror reflection. 

nomena at different frequencies. In figures |^a and |^b we depict the propagation of the 
shear elastic wave in the homogeneous plate and in the presence of the periodic structure, 
respectively. The spatial energy repartition and the cylindrical phase shift in the homo¬ 
geneous plate express the consequence of a point source generating shear vibrations. The 
non-dispersive feature of shear waves makes them travel with the same group and phase 
velocities and hence maintains the same behaviour regardless of frequency. The energy iso¬ 
surfaces allow us to follow the energy decay in the homogeneous plate. In £gure|^bl, in the 
presence of the periodic structure, the 28 kHz-frequency corresponds to the directional band 
gap for shear waves (second mode in figure |^c). This leads to a self collimation (or endo¬ 
scope effect). Through the phase shift analysis, we notice that the curvature of wavefronts 
is reduced thanks to the periodic structure and the energy isosurfaces are parallel over a 
longer distance. The energy repartition corroborates this observation. The most interesting 
effect can be seen at frequency 32 kHz. Indeed, this frequency represents the aforementioned 
inertial twisting mode (figure |^b). The sub-wavelength resonators convey collectively the 
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Figure 5: a) Transmission loss through a seven-column platonic crystal made of rotators with 
tilted ligaments. The solid blue curve describes the case when the rotators are placed in the right 
direction to convey the elastic waves and the dotted red when the rotators counteract its energy 
flow. The two curves fit perfectly except for frequencies close to resonance, b) Example of crystal 
response with counteraction effect at 32 kHz to be compared with Figure 4-b2. 

energy towards the point source image. At the crystal exit, we can clearly notice the change 
in wavefront curvature (from convex to concave) and the concentric energy isosurfaces cor¬ 
roborated by the magnitude cartography. Finally, we evaluated the well-known perfect 
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mirror effect caused by the complete band gap for shear elastic waves. 

In hgure [^a, we quantify the energy loss through the periodic structure collected from 
a sensor placed on the opposite side to the point source. The numerical simulations have 
been performed in both situations, when the point source is placed on the left side (red 
curve) and the right side (blue curve) of the periodic structure. One can hrst notice that 
the two curves retrace perfectly the position of the dispersion-based complete band gap. 
Furthermore, the red curve hts entirely the blue one except in the resonance region between 
30 and 40 kHz. We were particularly interested in the energy steering at the resonance 
frequency. In contrast to hgure |^b2, the rotators’ movements counteract the energy how 
and destroy completely the wave’s fringes (hgure 5-b). 

In conclusion, we analysed the propagation of elastic waves in a platonic crystal with 
inertial resonators. It is demonstrated that such a structured plate exhibits large band 
gaps due to the presence of heavy resonators linked to the host plate by tiny ligaments. 
In addition, the hat band representing the resonance frequency can be used to achieve new 
functionalities such as focusing. To do so, the ligaments should be tilted in order to facilitate 
the energy transfer of the twisting mode to the continuum, and vice versa. Since the platonic 
crystal is of constant thickness, the overall device could be implemented easily and might 
oher an alternative to free standing plates. Micro-Electro-Mechanical System applications 
could be targeted. 
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A. D. are thankful for the ERG funding through ANAMORPHISM Starting Grant. 
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